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1 Introduction 



The existence of a gluonium, which is a color singlet of the bound state of two or more 
gluons, is predicted by all the models of quantum chromodynamics (QCD), the lattice 
models |1] , the sum rules Q , the bag model and the effective Lagrangian approach || . 
Experimental searches of these states are in progress, see, e. g., for the reviews in ref. ||. 
It follows from the analysis of obtained results that the most likely candidates for glueballs 
are the following meson resonances ||: 



Table. The eventual candidates for glueball. 
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In connection with that it is very important to describe the gluonium spectra theo- 
retically. Attempts have been made earlier to consider gluonium in the framework of the 
potential model with massive structure gluons || [L5j analogous to the non-relativistic 
description of the quark- ant iquark system, ref. [|T6|j . 

At present, the relativistic single-time approach [pH-[]19[] is used widely for the descrip- 
tion of two-particle systems like quarkonium. The necessity of allowance for relativistic 
effects is caused by the fact that in many cases the contribution of the relativistic cor- 
rections is of the same order as the contribution of the non-relativistic Hamiltonian. In 
the present work, the quasipotential equal-time approach is employed for the description 
of two-gluon bound states [] consisting of the structural gluons which are described by 
the six-component Weinberg's wave functions pi|-P4|. The results of the papers []25], p6| , 
devoted to the covariant three-dimensional description of the composite system formed 
by two particles with the S = 1 spin, are used. The remarkable feature of our formalism 
is the locality of the corresponding quasipotential in the Lobachevsky momentum space. 
This is achieved by the separation of the kinematical Wigner rotations and "resetting 



1; 



all spin indices to the one momentum, for details see refs. Jl8|, ^7], [28|. Moreover, the 
quasipotential for interaction of two vector particles is the same as the quasipotential for 
interaction of two spinor particles with corresponding substitutions accounting for the 
spin difference and the normalization. 



lr The two-gluon bound system has recently been described on the basis of the Bethe-Salpeter equation, 
ref. |20|. 



1 



The quasipotential equation with the one-boson exchange potential, obtained in Sec- 
tion II, is reduced to the finite-difference partial-wave equations in the relativistic config- 
urational representation (RCR) in Section III. As was shown earlier ||29|| , it is possible to 
develop the relativistic analog of the WKB methods in the RCR. This method has been 
successfully used to find the quarkonium mass spectrum The use of the quasiclassi- 



cal quantization condition for relativistic two-particle states |T^, a], see Section IV, also 
allows calculation the gluonium energy levels. 



2 Spin structure of the relativistic potential for two- 
gluon interaction in the momentum representation 

At present, the gluon could be described as a massive particle with dynamical mass 
appearing due to the existence of color charge and self-interaction. This fact permits one 
to eliminate some contradictions in the results of calculations of the proton formfactor 
and the effective coupling constant as{q 2 ) on the basis of QCD (see in this connection 
ref. 0). 

Therefore, we begin by considering the quasipotential for two-gluon interaction in 
the momentum representation as that of gluonium consisting of the structure massive 
gluons with the intermediate interaction of the gauge massless gluon. The corresponding 
diagram describing this process is drawn at Fig. 1 . The Feynman matrix element 
< P11P2] o"2 I I ki, k 2 ; vi, v 2 > corresponding to this diagram is considered to be 
the quasipotential, V {2) = f {2 \ 




P2i a 2 k 2 , v 2 

To find the form of the quasipotential it is necessary to know the Feynman rules for 
the vertices of interaction of the structure gluons and the massless gauge gluon. 



In ref. [21] the attractive 2(25' + 1)- formalism for the description of particles of spin 



5 = 1 has been proposed. As opposed to the Proca functions, which transform according 
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to the (~, 5) representation of the Lorentz group in the case of S — 1, the spinor functions 
are constructed via the representation (S, 0)©(0, S) in the Joos- Weinberg formalism. This 
way of description of higher spin particles is on an equal footing to the Dirac description 
of spinor particles whose wave functions transform according to the (^,0) © (0, \) repre- 
sentation. The 2(2S + 1)- component analogues of the Dirac functions in the momentum 
space areQ 

U(v)- — ( DS ^^ \ (2 1) 

for the positive-energy states; and 

V(v)- — ( D*(a(p)C-i)C \ 

for the negative-energy states, with the following notations: 

p + M + (ap) 

«(PJ = / =» c = ( 2 - 3 ) 

pM (po + M) 

and D s (y4) = D^ S '°\A) is the Lorentz group representation by matrices with [2S + 1) 
rows and columns^. In the case of S — 1, one has 

(aW) = x + if + (, 4) 
In spite of some antiquity of this formalism, in our opinion, it does not deserve to be 



retired. Recently, some attention has been paid to this formalism p2|-p6 



In the articles b,h-j] the Feynman diagram technique is discussed for the vector 
particles in the above-mentioned six-component formalism for quantum electrodynamics 
(QED). The following Lagrangian: 

C QED = ^(x)r^V M V^(x) - M 2 ^(x)^(x) - \f^ u + 

+ ^F^(x) l5 ^(x) + ^^d a F^(x)^ >a pV^(x) (2.5) 

has been used there. In the above formula we have = —id^ — eA^; F^ u = d^A u — 
d v A^ is the electromagnetic field tensor; A^ is the 4- vector of the electromagnetic field; 
ty, \I/ are the six-component wave functions (WF) of the massive vector particle. The 
following expression has been obtained for the interaction vertex of a vector particle with 
a photon Pj], m : 

- eT af3 (p + k)p - ^l5,apqf3 + -^pl^ap^qpq^P + k)„, (2.6) 

where T a p = + 5 a/ 3; 7 a/ 3; 75,0/3; "Yq^p^v are 6 © 6- matrices which have been considered 
in ref. [fH|, b,g] 

2 These functions obeys the orthonormalization equations, C/ + ' (33)744 f (p) — 1 and analogous equation 
exists for V(p), the functions of negative-energy states. 

3 The technique of construction of D S (A) could be found in [|TJ k]. 
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5ij S{Sj SjSi 



djj SiSj SjSi 



7i4 = 74i 



iSi 
-iSi 



744 



1 

1 07' 



(here, Si are the spin matrices for a vector particle), 

l5,af3 = ifr/annPft], (2-7) 
l6,aP,fj,u = [japn/3v]+ + 2SaiM8pv — [yi3ft,lav]+ — 28i3 ft d a u, (2.8) 

e is electron charge, A and k are the quantities which correspond to the magnetic dipole 
moment and the electric quadrupole moment, respectively; M is the vector particle mass. 

In the case of interaction of two structure gluons we also use the 2(25 + 1)- component 
formalism. Since the gluon 2(2S + 1)- dimensional WF can not be directly introduced 
into the Lagrangian by means of the standard procedure of lengthening the derivative 
(covariantization), we add the terms defining the structure gluons manually: 

C^g D = gV^ v V v g - M 2 g gg (2.9) 
into the commonly used QCD Lagrangian: 

C QCD = iffT" V^g - m q qq - \g%G^ . (2.10) 

Where 

V ll g(x) = d ll g-ifT^B^g (2.11) 

and g, g are the 2(2/9+1)- dimensional WF of the color octet, T° d - are the SU (3) generators 
in the adjoin representation. Employing the technique of functional integration and using 
a scheme which is analogous to the qqg- vertex case, we suggest that the interaction vertex 
of two structure gluons with the gauge gluon has the following analytical formQ: 

if-rip + k)v{ThU (2-12) 
without taking into account the multipole momenta (compare with the QED expression 



Using results of refs. f25|, |28|j, let us represent the Feynman matrix element correspond- 
ing to the diagram of one-gluon exchange (see Fig. 1) as: 

< pi,p 2 ] o- 1 ,a 2 \V (2) \k 1 , k 2 ; vi, v 2 >=< Pi,p 2 ; a 1 ,a 2 \T ( - 2 ' ) \k 1 .k2; v x , v 2 >= 



&ip i Vip i^ih — 1 

^£3 [V- 1 ^ k 2 )}D^ {V-\A V , k 2 )} , (2.13) 



4 We did not allow for the term S^ v (p + fe)„ which is the result of the auxiliary condition (Klein- 
Gordon equation) to the Weinberg equation. In the Kadyshevsky's approach all the particles, even in 
the intermediate states, are on the mass shell. 



4 



where 



After some calculation, using the formulae of ref. |2j 

UM = S p 4J a (0), S^S n = S S( _x (v-\a p , k) 

S^^l^uPuSp = 744 {Pn ~ JsWuip)) , 

WM ■ d (v-\A f , k)) = d (y-^Ae, k)) ■ U,(k) + J^+ M) p v w v {k) 
k,w,(p) ■ d (y-\Ap, k)) = -d (y-\k p ^ k)) ■ p.W^k), 

we come to the 4- current of vector particle^: 

f^(p, k) = -f^ p {(p + k), + ±W,(p)(SA) - -L(SA)W,(p)}z Up 



and to the quasipotential 

V (2 \kHP,P) = -3f 2 { 
10! + S 2 ) \pA 



p Q (A + M) + (pA)} -M 3 (A + M) 



M 3 (A — M) 



+ 



p (A + M) +pA 
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1 St \pA\ S 2 pA 



M 3 
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As used in the earlier works |18|, |28|], we have 

A = A;Jk = k(-)p = k-^(k - 



p kp 



Po + M' 



A = (A^) = (k oPo - kp)/M } 



(2.14) 

(2.15) 
(2.16) 
(2.17) 
(2.18) 

(2.19) 



(2.20) 



(2.21) 



(2.22) 



where p, k are the covariant generalizations^ of the vectors of particle momenta in c.m.s., 
before p\ = —p 2 = p and after k\ = —k 2 = k interaction; £*, £ are the analogues of Pauli 
spinors and D^g is the Wigner matrix of the irreducible representation of the rotation 
group, which has dimension equal to (25 + 1) with the following form: 



D^{v-^k)} 



2M(p + M)(k + M)(A + M) 



pk 



+ 



+ 



(p + M)(k + M) - kp 



2? 



(p + M)(k + M) — kp S pk 



2{S[pk}} 2 } 



(2.23) 



^Wfj_ is the Pauli-Lyuban'sky 4- vector of relativistic spin. 

3 We omit the circles above the covariant generalizations of the momenta, as opposed to 
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The expression ( |2.20| ) shows the advantages of the 2(25+ 1)- formalism, since it looks 
like the quasipotential for the interaction of two spinor particles with the substitutions 
75— rk — \ =>- ~ko an d S =>- a. 

2m(A-m) A 2 



3 System of two-particle partial-wave equations in 
the relativistic configurational representation 

The transformations into the relativistic configurational representation (RCR) have 
the following form: 

V(r,n;p) = — L- / dtl A £*(A; f)V(A,p), (3.1) 



for the quasipotential and 



(2tt) 3 
1 



(3.2) 



(2tt) 3 

for the WF. The integration measure is 

dQ p = d 3 p/^fl+p 2 /M 2 (3.3) 

It is the invariant measure on the hyperboloid, Pq — p 2 = M 2 . The system of functions 
n, r) is the complete orthogonal system of functions in the Lobachevsky space, 

f Po -pn w _ irM 



M 



(3.4) 



The physical meaning of the parameter r is discussed in details in ref. [[31 



As a result of earring out this transformation to the RCR we arrive at the following 
quasipotential: 

V(r,n;p ,p) = Vi(r,p Q ) + V 2 (r,n;p ,p), 



(3.5) 



with 



3PI 



1 



V^po) = -3f \ (8p 2 - AM 2 )V Yuk (r) + - l)-5(r 2 + —)+ 



— — 5(r H ) + 

M 2 r y t m 2; m 2 



M 2 



(Si £2 



2p 2 11 I 11 4 . 



(3.6) 



and 



V 2 (r,n;p 0} p) = -3/ 



_3f2 / A(^(^ S ( r ) _ 5l2jB ( r ) 



IM 2 



+ 
+ 
+ 



M 2 
2ip 



{S 1 L)(S 2 L) + {S 2 L)(S 1 L) 



:B(r) 



M 2 



(pn) 2 B(r) + 



(pn) 



4rA(r) - — C(r) 



M 

(jm)-B(r) — ^j^~C(r 



(SL) 



4 P° A( \ _i 



M 2 



M 3 



(3.7) 



where S = Si + S 2 , L = [p x f], S 12 = 3(S 1 n)(S 2 n) - (5152). 
Next, 

V Yu k{r) = -^-cth{rMTT), (3i 



and 



A{T) = r(r4/M)) ' Wr> - <3 ' 9) 

B(r) = (r + (i/M))U(2*/M)) ^' (3 - 10) 

_ 1 (r-(,:/Af)) l 4 
° W " (i/M) r + (i/M) r [r + hP> 4r (iU > 

From the above equations we can see that the quasipotential is separated into two 
parts, Vi(r;p ), which does not depend on the direction of the " relativistic coordinate" 
vector, and V 2 (r, n;p ,p), which depends on n over the structure (pn), \p x n], S\ 2 • 

After the transformation of the quasipotential equations: 

(M - 2 Po )V aia2 (p) = (2vr)- 3 J2 I dn k v:^ 2 (k,p)^ ulU2 (k), (3.12) 

U\U2 

in the RCR by means of the formulae ( |3.1| , |3.2|) , we can be convinced that Vi(r; p ) describes 
the local type of interaction and V 2 (r, n;p ,p) enters into equation (|3.13|) by the non-local 
way. 



x ^(p;^,^)^(p;^i,^i)V2^ 2 (ri,niAp;po,^^^ 1 ^(^i) ) ( 3 - 13 ) 



where the unit vector is [|32l 



n Ap = [Mn - p(l - -J^—)]/(p - pn). (3.14) 
Po + M 

It is still possible to localize the spin-orbit part, some terms of the tensor interaction and 
some other terms entering in V 2 (ri, Uia ;po,p). 
We have the following structure in Eq. (|3.13| ): 

(p-n 1Ap ) = M 2 /(p -prii) - Po , (3.15) 

[p x niAp] = M[px ni] /(p - pHi), (3.16) 



and 
where 



(5 ? ini Ap )(S ? 2 n 1A p) = Z\ + Z%, (3.17) 

Z T X = M 2 (Sini)(S 2 fii)/ (p — pni), (3.18) 

M 2 1 -» -. 

# = (ffQ _^ )8 {-jfW(^ + 0M(3»*i)] x 

x ( 1 -^^) + 4(^(^( 1 -^ i 7) 2 }- (3-19) 
Po + M p + M 
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The localization procedure for ( |3.15|jrT6"|) and the first term of Z\ is produced by 
means of the following equation: 



M OTi Pq — prii 



As a result we obtain 



(3.20) 



(3.21) 



The quasipotential of Eq. (|3.13|) is presented as a sum of six component: 
V(f;p 0l p) = V C + V L s(LS) + V s 0xS 2 ) + VtS 12 + V LL L 12 + V Sp P 12 , 

where 

L 12 = ^{(S ? 1 L)(S ? 2 L) + (S ? 2 L)(S ? 1 L)}, 
P12 = (SiP)(S 2 p); 



(3.22) 

(3.23) 
(3.24) 



and 



LS 



+ 



Vo 



Vll 
Vs P 



3/ 2 {(8p - 4M*)V Yuk (r) + (M _ 1) V + 4) + 4 V + T77 



+ 2 



p 2, 
M 2 



2 (9 tP" 



4rA(r) - ^C(r) 



5(r) + 



M 2 
-6 

2zpo 
M 



r 



M 2 M 2 r 



M 2 ' 



tzll! exp( _^ ) _2 | ,(r-i) !x 



(r - -M 2 



exp( 



Mdr' M r 2 
i d , 



r 2 ""^ M9r j Po 



x 



M (r + (i/M)) 5 



M 



-3/ 2 

M 2 
-3/ 2 

1 1 S , 2 1 x 1 1 r , 2 1 



+ 



r(r-{i/M)) id 

(r - (2i/M)f . 2% d . 2f 



Mdr' M 2 



r + (i/M)) 2 
B(r 



B(r) 



-3/ 2 
-3/ 2 



(r-(2i/M)) 2 2* 9 

5 ex P(-i7^r) 5 ( r ) 



M<9r' 



M 2 (r + (i/M))(r + (2i/M)) 



-B(r) 



Zf{ W2 B{r)}. 



(3.25) 



(3.26) 



(3.27) 
(3.28) 

(3.29) 
(3.30) 
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It is noted that we have assumed the corresponding finite-difference operators in po 
and p for these calculations. However, the use of eigenvalues instead of these operators is 
a good approximation, which does not reduce to the non-relativistic models. 

After carrying out the partial-wave expansion for the WF 

* {S) (r; °) = ^J2 RMr) , (3.31) 

JIM 

the three-dimensional quasipotential equation is rewritten in the system of radial equa- 
tions for S — 0, S — 1 and 5 = 2 

(M - 2H e )Rjes(r) = £ % J s/ , s ,{r; P(h p)R M , s ,{r). (3.32) 

VS' 

Here, 

H e = Mcoshi ——) + -sinhi——) + \ r J exp — — , 3.33 

y Mdr> r y Mdr' 2Mr 2 FK Mdr h v ; 

and ^ 

VeS'M 1 "' Po>&) = J du H njf M (n)V(f; p , p)ttf^ M {n). (3.34) 

One can check that due to the complete accordance of the relativistic spin structure of 
the quasipotential ( |3.22| ) to the spin structures which are used in non-relativistic models, 
the matrix elements ( |3.34j ) have the same form as when non-relativistic interaction of two 
vector particles is considered. The only differences are in the internal expressions of Vc, 
V T , Vsl, V s , Vll and V Sp . 

Thus, we have for a singlet state 

(M - 2H e=J )R e=J (r) = (V c - 2V S - ~ J( J + l)V LI )j Re=j(r); (3.35) 
for a triplet state 

(M - 2H e=J+1 )R e=J+1 (r) = (v c - ( J + 2)V SL -V s - ^j^V T + 



J + 2- \ 3JJ(J + 1), 

+^~Vll) Ri=j + i(r) + 2J+1 VrRt=j-i(r), (3.36) 

(M - 2H e=J ) R e=J (r) = (v c - Vls ~ % + Vr + [~ - J {J + 1)]Vll) Re=j(r) } 

(3.37) 

(M - 2H e=J ^)R e=J ^(r) = (v c + (J~ 1)Vls ^f^Vr- 

J-l- \ 3 Jj( J + 1) , 

Vll )Ri=j-i(r) + _ T , ,, V T R £=J+1 (r); (3.38) 



2 J y ' 2J + 1 

for 5-plet state 

(M - 2H e=J+2 )R e=J+2 (r) = (v c - 2(J + 3)V LS + V S - ^ + jj V T + 
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+{J + 3) 2 V LL )R e=J+2 {r) + 



6J(J + 2) J(2J- 1)(2J + 1) 



2J+ 1 

(A* - 2H £=J+1 )R e=J+1 (r) =[V C -(J + A)V LS + V s 



2J + 3 
J-4 



2J + 1 



3J + 8, 



3^(J-l)(J + 2) 
2J+ 1 



+ 



/6J(J + 2) J (2 J - 1)(2J+1) A 
(.M - 2H e=J )R e=J (r) = v o ; ■ - ;VV oT V o -V» =J+2 (r) + 



(3.39) 



(3.40) 



2J+ 1 



2J + 3 



\ C LS S (2J-l)(2J + 3) T 



5 1 
2 ~ 3 



j(j + r 



Vr r X 



/6(J-l)(J+l) v /(2J+l)(2J+3) 
xi? £=J (r) + J! — — * — — V T R e=J ^(r) 



2J+ 1 



2,1 - 1 



(3.41) 



(M - 2H e= j_ 1 )R e=J ^(r) = [Vc + ( J - 3)\> L5 + K< 



J + 5 
2J + 1 



TV 



3J-5 T> \ _ . . 3^/(J-l)(J + 2) 

-Vll ) Re=j-i{r) + - 1 — 77T —^ V T R e=J+1 {r), 



2J+1 

(M - 2H t= j- 2 )R^j-2(r) = (% + 2(J - 2)T> XS + V 5 - ^jlf ^+ 



(3.42) 



, , ,/6(J-l)(J + l) J(2J + l)(2J + 3) 
-(J - 2)V LL ) i?, =J _ 2 (r) + V , T A - V . T \ -V T Ii, ,(/■) 



2J+ 1 



2J- 1 



(3.43) 



The two-gluon bound states are to have the positive C- parity. Therefore, we have 
restrictions on J in ( |3.35|4T43l ). 

Finally, let us remark that we have neglected the last term of ( |3.22| ) following the 
authors of ref. [091. 



4 Quasiclassical condition for a quantization in RCR 

The non-tied partial-wave equations can be re-written in the following form: 

-8 i A -8 X 2 f(f + 1) -8 

[cosh(i\—) + -srnh{f\-) + L± I exp(.A-) - X{r)]R,{r) = 0, (4.1) 

W - V(r) 

X( r ) = (4.2) 

V ; 2M ' V ; 

where W = 2M + is the mass of bound state, M is the gluon mass, is the binding 

energy. We consider V(r) as a sum of the potential V con f describing the confinement of 

a gluon in a meson and the corresponding matrix elements of the quasipotential in the 

system of the Eqs. ( |3.32| ). In the expressions ( |3.25| - |3.30| ) it is possible to neglect the 
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small image additions which are proportional to i/M(= ih/Mc), because the bound state 
spectrum usually forms on distances of the order r » X = H/Mc. 

The quasiclassical condition for a quantization of the two-particle relativistic systems 
has the following form : 

/ r+ dr'ArchX k (r') = Xrr(n + -), (4.3) 

where 

X A (r)=X(r)[l + (A^) 2 ]- 1 / 2 , A = i + 1, (4.4) 

and the integration limits are determined from the equation 

X A (r ± ) = 1. (4.5) 

If restrict ourselves to the cases of the simplest potentials, namely, V(r) = or and 
V(r) = or 2 , the quantization condition is shown in ref. |34j to take the form of Q4.6| ) and 
(|4.7|) ) correspondingly, 

a £3 

Y cosh v — sinh y = -A7r(rH 1 — ), (4.6) 

2Mc 2 2 4 ; v ; 

i — ^ ^ 

2 v &^TT[K(tan/ iX /2) - S(tan/ix/2)] = \/^^M n +2 + 4 ); (4 ' 7) 

where cos/i x — M^/2M. Following this technique, the energy levels of gluonium states can 
be obtained. The numerical results for the case of the confinement potential are presented 
in ref. [26|. Investigation of the influence of the Coulomb term as well as the relativistic 
corrections are in progress. 



5 Conclusions 

In the present work the formalism for consideration of the two vector particles, gluons, 
has been constructed on the basis of the 2(2S + 1)- dimensional description of the WF. 
The form of the relativistic two-particle single-time quasipotential equation was found. It 
turned out that the quasipotential of this equation coincided with the quasipotential for 
interaction of two spinor particles in the second order of perturbation theory. This fact 
shows the advantages of the 2(25' + 1)- component formalism. 

It is possible to employ the relativistic generalization of the WKB method to find 
the gluonium spectrum in the case when the total spin of the system is equal to zero. 
In following works we will use the obtained system of the partial-wave equations for 
investigation of the contributions of the spin-spin and spin-orbit interactions to the energy 
of the gluonium states. 
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